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Model-Independent Semileptonic Form Factors Using Dispersion Relations 

C. Glenn Boyd, Benjamm Grinstein and Richard F. Lebed 
Department of Physics, University of California, San Diego 
La JoUa, California 92093-0319 

We present a method for parametrizing heavy meson semileptonic form factors 
using dispersion relations, and from it produce a two-parameter description of 
the B —>■ B elastic form factor. We use heavy quark symmetry to relate this 
function to B D*W form factors, and extract |Vcb| = 0.0355^q;qq25 from 
experimental data with a least squares fit. Our method eliminates model- 
dependent uncertainties inherent in choosing a parametrization for the ex- 
trapolation of the differential decay rate to threshold. 



1. Introduction 

A non-perturbative, model-independent description of QCD form factors is a desirable 
ingredient for the extraction of Cabibbo-Kobayashi-Maskawa parameters from exclusive 
meson decays. Progress towards this goal has been realized by the development of heavy 
quark symmetry[jl|, which relates and normalizes the B —>■ D*li/ and B —>■ DW form factors 
in the context of a -^^ expansion, where M is the heavy quark mass. Previous talks 0] have 
described how this normalization is used[|^-|^ to extract the value of the CKM parameter 
\Vcb\ by extrapolating the measured form factor to zero recoil, where the normalization is 
predicted. 

This form factor extrapolation, necessary because the rate vanishes at zero recoil, 
introduces an uncertainty in the value of \Vcb\ due to the choice of parametrization. Esti- 
mates of this uncertainty obtained by varying parametrizations suffer the same ambiguity. 
This ambiguity could be eliminated if one had a non-perturbative, model-independent 
characterization of the form factor in terms of a small number of parameters. 

In this talk we describe recent work[^|-|^ in which we use dispersion relations to derive 
such a characterization and apply it towards the extraction of \Vcb\- Our characterization 
of the B elastic 6-number form factor uses two parameters and has 1% accuracy over the 
entire physical range relevant to the extraction of \Vcb\- What is important is not that we 
effectively find a quadratic parametrization of the Isgur-Wise function, but rather that we 
have determined the associated uncertainty. 
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The outline of the talk is as follows. We first review the expected achievable accuracy 
in the extraction of \Vcb\ from semi-inclusive semi-leptonic decays in Sec. 2. In Sec. 3 
we review a well-known method for using QCD dispersion relations and analyticity to 
place constraints on hadronic form factors, and show that a parametrization of the B 
elastic form factor F[q^) in terms of two parameters is accurate to 1% over the relevant 
kinematic region. In Sec. 4 we use heavy quark symmetry to relate F{q^) to the Isgur-Wise 
function, which describes the form factors for B D*W in the infinite quark mass limit. 
We make a least squares fit to CLEO0, ARGUS [g, and ALEPHg data using \Vcb\ and 
our two basis function parameters as variables, and present our results. Reliability of the 
method is discussed in Sec. 5, implications for \Vub\ and B K*^ are discussed in Sec. 6, 
and a summary in Sec. 7. 



2. Digression: \Vcb\ from inclusive semi-leptonic decays. 

Several competing methods for the determination of \Vcb\ were described in previous 
talks0. We would like to discuss briefiy the theoretical limitations in the determination 
of \Vcb\ from inclusive semi-leptonic decays. The interpretation of the measurement of the 
inclusive semileptonic rate r{B — s> Xcii^i) relies on our ability to calculate the rate from 
first principles. Using a Heavy Quark expansion one can show two things 

(i) The leading term (in l/nib) is given by the parton decay rate r(6 ciui) 

(ii) There are no first order (in 1/mb) corrections to the previous statement. 
In a l/nib expansion one may write p!0| , p | 



(2.1) T{B ^ X,^Ul) = To[A{x)r^ + BK{x)Kb + BG{x)Gb + 0{l/ml)]. 

There is a lot of notation to explain here. First, TqA{x) is just the parton decay rate, with 
To = G'|.m^|Vc&P/1927r^, x = ml/ml and A{x) = 1— 8x-t-8x^— x* — 12x^ Inx is a kinematic 
factor from the phase space integral, r] is a correction factor |Tn] from perturbative QCD, 



(2.2) ^ = 1 - O ^^'^(^) - (f^oAf^lix) + ■■■) + 0{a'). 

The one loop function A^^\x) can be computed analytically and equals 1.7 at a; = (0.3)^ 
(2.8 at X = 0). The full two-loop computation is not available, but the part proportional 
to /3o, the one-loop beta function, has been computed: A^^-* = 1.7 at x = (0.3)^ (3.2 at 
X = 0). 
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The next two terms are the 1 / corrections, with 
(2.3) Kh = ^{B\b^^b\B) and = ^{B\b^^^^^^b\B) 



and the corresponding kinematic factors Bk{x) = —A{x) and Bq[x) = 3 — 8a; + 2Ax'^ — 
24:X^ + 5a;^ + 12a;^ In a;. 

There are three main sources of thoretical uncertainties: 
• Quark masses enter the rate in the combination m^A(x). This is less sensitive to 
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uncertainties in than m^ if mc is fixed by 



(2.4) ms - rriD = mt - mc + mb{Kb + Gb) - rndKc + Gc) 

For a fixed guess oiKb,c and Gh,c Wisefig (Shifman|ll|]) finds that /\\Vcb\/\Vcb\ = 10% 
(3.2%) for Amb = 0.5 GeV (0.2 GeV). 

Perturbative QCD corrections are stiU not fuUy computed at two loops. This is not 
a limitation in principle, except in the case that the expansion does not converge 
quickly.! With x = (0.3)^, /3o = 9 and = 0.2 the three terms in (2.2) are0 



T] = 1 — 0.11 — 0.06 (77 = 1 — 0.15 — 0.11 for x = 0), which appears only marginally 
convergent. 

Non-perturbative corrections come in through the parameters K and G. The latter 



is well know as it determines the B* — B mass difference. Shifman ||ll|| estimates an 
uncertainty of A\Vcb\/\Vcb\ = 5.6% from Kb = 0.024 ± 0.009. 

I close with a digression on a method related but not the same as the heavy 
quark expansion. One may argue that r{B Xcii^e) = r(^ ciui) in the limit 
mb — m,c ^ Aqcd ^ ^c,f)- The argument is straightforward. In this limit only the 
D and D* resonances are kinematically allowed, so they saturate the semi-inclusive rate. 
The rates into these resonances are determined effectively by the Isgur-Wise function at 
V ■ v' = 1, where it is normalized to unity. The resulting rate is precisely r(6 ciui), 
provided one does not differentiate between quark and meson masses. 

In the limit, the rate for r{B — ^ Xiui) ~ (m^ — mx)^, for X = D, D*. This can be 
used to 'explain' the lower sensitivity of r{B Xcii^i) to mb when mc is fixed in terms 



^ The expansion is asymptotic, so it does not converge. The asymptotic expansion has an 
accuracy equal to the size of the term at which the expansion is truncated. 
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of mi, using the known value of — mD . But the real question is whether this could be 
the basis for a systematic expansion. Note, however, that 



(niB - mnf = (m^ - m^f{l + K + Gf + ■ ■ ■ (2.5a) 

(niB - niD*)'' = {rrib - m,)^(l + K - G)'' ( 1 + IOG ""^^""" ^ + ■ ■ ■ (2.56) 

V mb-rricj 

These correction factors will enter the relation between inclusive and partonic widths. The 
expansion parameter in Eq. (p.55|) is poor: Qin±±lP- 



mb—rric mf,—mc mj, 



3. The Analyticity Constraints^ 

Of primary interest is the form factor F, defined by 

(3.1) {B{p')\V^\B{p)) = Fiq^)ip+p')^, 

where = by^b, and = {p — p'Y- Crossing symmetry states that the form factor F is 
a function of a complex variable, q^, which gives the matrix element in Eq. (3.1) for real 
negative q^ , and the matrix element {Q\V^\BB) for real q^ > 4m^. 
Consider the two-point function 

(3.2) i J rfW^^(Ty^(a:)yJ(0)) = {q^q, - q^g^^Mq^) 

In QCD the structure functions satisfy a once-subtracted dispersion relation: 



(3.3) 



dq'^ 



1 ['^ , imn(t) 

' dt ^ ' 



'2 = -Q2 



Jo {t + Q 



2^2 ■ 



The absorptive part Imn(q'^) is obtained by inserting real states between the two currents 
on the right-hand side of Eq. (3.2). This is a sum of positive definite terms, so one can 
obtain strict inequalities by concentrating on the term with intermediate states of BB pairs. 
For far from the resonance region (Q^ -|- 4ml ^ "^b^QCo) the two-point function can 
be computed reliably from perturbative QCD. We set = which for large b quark mass 
is far from resonances. One thus obtains an inequality of the form 

POO 

(3.4) / dtk{t)\F\'^ <1, 

Jam% 



This section based on refs. |^,12-14] 
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where the function k{t) is the ratio of the kinematic factor on the right hand side of 
Eq. (3.3) to the QCD calculation of the left hand side. 

A key ingredient in this approach is the transformation that maps the complex 
plane onto the unit disc \z\ < 1: 



(3.5) 



1 + z 



AMI 1 



In terms of the angular variable e*^ = 2, the once-subtracted QCD dispersion relation may 
be written as 

(3.6) ^ r d9me'')F{e^')\'< - , 

where the weig hing function <P{e'^) = k{t{e))dt/d6^: 



(3.7) <Piz) = ^^^'-^{l + zrVT^z. 

Here n/ is the number of light flavors for which SU{nf) flavor symmetry is valid; we take 
Uf = 2. Perturbative corrections to the dispersion relation are incorporated in p, which 
has been computed [p!5[] to 0{as), p = 1 + 0.73as{mb) ~ 1.20. 

Note that (j){z) is analytic in |z| < 1, while poles of F inside the unit disc originate from 
resonances below threshold and cannot be ignored[0. In this talk we only consider the 
effects of the resonances Ti 2,3. Although cuts below threshold should be considered, they 
are expected to have smaller physical effects. A simple but effective trick[j4] eliminates 



the poles with no reference to the size of their residues but rather only their positions {i.e. 
masses). The function 

(3.8) P{z)- i^-^i)(^-^2){z-zs) 



(1 - Ziz){l - Z2Z){\ - Zsz) 

where the Zi correspond to the values q^ = M^-. , is analytic in 1^1 < 1 and satisfies 
|P(;2)| = 1 for \z\ = 1. The function P{z)(p{z)F{z) is analytic on the unit disk, and obeys 

(3.9) ^[ de\P{e''^)(P{e''^)F{e'^)\^ 



2tt Jq n 
It follows that the QCD form factor may therefore be written as 

(3.10a, b) F{z) = , V a^^" with V |a„|2 < - 

Pizjmz) ^-^ ^-^ TT 
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Since S-number is conserved, F{0) = 1, so that ao = P(0)^(0). 

In the next section we wiU use heavy quark symmetries to relate F to the form factors 
for B —>■ D*li/, where the physical kinematic range is < z < 0.056. Therefore, the form 
factor in (3.10a) converges quickly over the physical region. Note that for this it is crucial 
that the coefficients be bounded as in (3.106). Retaining only ai and 02 in Eq. (3.10a), 
we have 

(3.11) F{z) = -^^-L^[P(0)(/)(0) + aiz + a2Z^] 

with a maximum relative error of ~ v^T7^(0.056)7P(0)(/)(0) ^ 0.01. 



4. Extraction of \Vcb\ 

4-1. Heavy Quark Symmetry Relations 

In the infinite b and c quark mass limit all the form factors for B Dlu and B —>■ D*li/ 
are given by one universal "Isgur-Wise" function. This allows us to apply the constraint 
on F to the particular combination of form factors actually measured, rather than deriving 
constraints for each form factor separately. For large 6-mass, the Isgur-Wise function is 
given by the form factor F, F{uj) = ^(cu), where uj = v ■ v' = pb ■ Pd/i^bi^d- There are 



short-distance matching and running corrections! 16 1 to the relation between ^(u;) and the 
B D* form factors. For example, for the vector current form factor g one may write 
g{z) = rjDF{z). This relation generally holds to order 1/M, but at threshold it holds to 
order l/M^I^. We take rjn = 0.985. 

4-2. Maximum Likelihood Fit 

Once the essential physics of QCD is incorporated into the calculation via Eq. (3.11), 
the maximum likelihood fit is simply an ordinary chi-squared minimization with parameters 
|Vcf,|, ai and a2. We normalize input data to a S lifetime [[T8[| of tb = 1.61 ps. Also, we 



rescale the ARGUS data to refiect a revised branching fraction for — > K~'k^. 

Table 1 shows the central values and 68% confidence levels for IK&I, ^i, and a2 from 
the various experiments. That a2 saturates the bound |a2| < 0.55 is not significant because 
its variance is large. Figure 1 shows the product of the best fit form factors with |Vc&|, 
superimposed with experimental data. At 90% confidence level, ai and a2 are consistent 
with zero, suggesting the dispersion relation may be saturated entirely by higher states. 
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\Vcb\- 10^ 


ai 


a2 


Expt. 


09 9+4.58 
■^^•^-5.89 


'-'•'-"-'-0.07 
'-'■-'--'--0.03 
'-'•'-"-'-0.04 




CLEOi 
ARGUS @ 
ALEPHg 



Table 1. Fit values for \Vcb\, ai, and a2 from the various experiments. 

The errors on |Vc&| in Table 1 are statistical only; the treatment of systematic errors 
depends both on our parametrization and a detailed understanding of the experiment. The 
error implicit in the variation over choices of parametrization, however, is absent. ARGUS 
examined the effect of varying over four possible parametrizations, which induced a spread 
of 0.012 in \ Vcb\, and was the major impediment in using heavy quark symmetry to obtain 
a model-independent extraction. 
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1.4 



Figure 1. Best fit values for 
CLEO (solid line), ARGUS 
The data for each experimei 
ization used in the text, is 
to reflect a revised branchin 



5. Discussion of the Method 

Perturbative and non-perturbative corrections to the dispersion relation (3.9) enter 
our analysis by modifying the value of the parameter p in Eq. (3.7). The extraction of 
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\Vcb\ is rather insensitive to such corrections. A change of p by ±10% changes the central 
value of \Vcb\ by less than ±0.3%. 

The truncation of Eq. (3.10a) to the first terms introduces an error proportional to 
p(o)Vo) Sn=Af+i '^nZ'^ < 1%, & bouud Valid over the whole physical region < 2; < 0.056. 

The application of the B B dispersion relation to B D*W decays relies on heavy 
quark symmetry. This is potentially the largest source of error, of 0{jj). We estimate 
such corrections by making a 20% change in the ranges of ai and 02, resulting in a 2% 
shift in the central value of \Vcb\- One can derive analogous dispersion relations for each 
of the B — > D*W form factors, sidestepping the need for heavy quark symmetry (except 

Another heavy quark correction arises at 0(^2) in the normalization of the Isgur- 
Wise function at threshold. The normalization of the form factor g{uj = 1) has been 
estimated to be g{l) = 0.96^, g{l) = 0.89|2l|], and g{l) = 0.93[|2l. We have included a 



QCD correction of 0.985, so to good approximation, this simply rescales the values of |Vcb| 
in Table 1 by 

3(1) 

There are other errors in our extraction that are not purely theoretical. The most 
pressing of these involve the binning of the measured rate against cu, smearing of uj in- 
troduced by boosting from the lab to the center of mass frame, and correlation of errors. 
Randomly varying input values of w in our least squares fit of the CLEO data by ±0.05 
changes the central value of \Vcb\ by less than 1%. A more thorough extraction can be 
done by the experimental groups themselves, using our parametrization in their analysis. 



6. Other Applications. 

Our parametrization of form factors applies to other heavy hadron decays, including 
B —>■ nW, with minor modifications. In this case the range of the kinematic variable is 
larger, < z < 0.5, so more coefficients a„ are needed for comparable accuracy. We expect 
six to eight will be necessary for accuracy of a few percent over the entire kinematic 
range, depending on the form of the actual data. 

Heavy quark symmetry relates the form factors for B nlu and D nlu. The latter 
is readily measured over < z < 0.3, and our method then allows a reliable extrapolation to 
0.3 < z < 0.5. This opens up the possibility of reliably extracting \ Vub \ from a measurement 
of Br{B nW). 
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An extraction of a model-independent lower bound on \Vub\ should be possible since 
small values of \Vub\ tend to wash out the nontrivial z dependence, while the a^'s cannot 
compensate because they are bounded from above. 

Any model of hadronic form factors must predict a form factor that is consistent with 
our parametrization. This is a severe test to pass, and serves as an effective discriminator 
for mo dels 0. 

Our parametrization may also be useful in the analysis of B ^ K*'j, by relating its 
amplitude to the form factor for D K*eu extrapolated outside the physical region [p3| . 

7. Summary 

The extraction of the CKM mixing parameter |V"c6| involves several types of uncer- 
tainties. Typically, these uncertainties are classified as 

(7.1) \Vcb\ = V ± {stat} ± {syst} ± {life} ± {norm} ± {param} 

where stat and syst refer to statistical and systematic experimental uncertainties, life refers 
to uncertainties in the B lifetime, norm refers to uncertainty in the value of the form 
factor at threshold, and param refers to uncertainty in the extrapolation of the measured 
differential rate to threshold. 

Not only the central value, but also the statistical uncertainty depends on the 
parametrization. For example, linear fits to CLEO data yield substantially smaller sta- 
tistical uncertainties than quadratic fits. Typically, quoted values correspond to the 
parametrization yielding the smallest statistical uncertainty, in effect throwing some sta- 
tistical uncertainty into the parametrization uncertainty, which remains implicit. Clearly, 
this does not improve the accuracy with which we know |V"cf,|. 

We have essentially eliminated the uncertainty in the choice of parametrization. This 
was accomplished in four stages. First, we used QCD dispersion relations to constrain the 
B —> B elastic form factor. Second, we derived a general parametrization of the B elastic 
form factor that automatically satisfies the dispersion relation constraint. This expression 
involved an infinite number of parameters bounded by J2'^=o < I- Third, we 
used heavy quark symmetry to relate the B elastic form factor to B D*li/ form factors 
and fixed the normalization at threshold. Over the entire kinematic range relevant to 
B — s> D*W, we showed that neglecting all but the first two parameters ai, 02 resulted in 
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at most a 1% deviation in the predicted form factor. Finally, we made a least squares fit 
of the differential B — > D*W rate to \Vcb\, oi, and a^. 

The results of this fit improve on all previous extractions in one important way: The 
uncertainty due to the choice of parametrization is under control. Our statistical errors 
are larger than many quoted values. This does not reflect an inferiority of our method, 
but rather quantifies uncertainties that were previously left implicit. Our averaged value 
from CLEO, ARGUS, and ALEPH data is 

(7.2) \Vcb\ X 10^ = 35.5tli (stat). 

An estimation of systematic uncertainties requires a detailed knowledge of the experiments. 

Our parametrization may be useful for other processes, such as 5 — > ttIP and B — > K*^ 
as well. For B ttW we expect to be able to extract a model independent lower bound 
on \Vub\, and a good measure of \Vub\ by use of heavy quark symmetries to infer the form 
factor from D nW by extrapolation. In addition, precision tests of QCD-predicted 
form factors are now possible; these should be useful as checks of QCD models and lattice 
simulations. 
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